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ABSTRACT 

We examine the practice of deriving interstellar medium (ISM) abundances from low-resolution spec- 
troscopy of GRB afterglows. We argue that the multi-ion single-component curve-of-growth analysis 
technique systematically underestimates the column densities of the metal-line profiles commonly ob- 
served for GRB. This systematic underestimate is accentuated by the fact that many GRB line-profiles 
(e.g. GRB 050730, GRB 050820, GRB 051111) are comprised of 'clouds' with a bi-modal distribution 
of column density. Such line-profiles may be characteristic of a sightline which penetrates both a high 
density star-forming region and more distant, ambient ISM material. Our analysis suggests that the ma- 
jority of abundances reported in the literature are systematically underestimates and that the reported 
errors are frequently over-optimistic. Further, we demonstrate that one cannot even report relative abun- 
dances with confidence. The implications are profound for our current understanding on the metallicity, 
dust-to-gas ratio, and chemical abundances of the ISM in GRB host galaxies. For example, we argue 
that all but a few sightlines allow for the gas to have at least solar metallicity. Finally, we suggests new 
approaches for constraining the abundances. 



Subject headings: gamma rays: bursts 
1. INTRODUCTION 



ISM: abundances 



With the launch of the Swift satellite (Cehrels et al. 

2004) , the rate of GRB detections has increased by more 
than an order of magnitude (Gehrels 2006). In turn, one 
has seen roughly the same increase in the detection of 
bright afterglows. Though a few high-resolution spectra 
were obtained pre-Swift (Castro et al. 2003; Fiore et al. 

2005) , Swift has enabled the nearly routine acquisition of 
high signal-to- noise ratio (SNR), high-resolution spectra 
(Chen et al. 2005; Prochaska, Chen & Bloom 2006). In 
addition, there are now over a dozen low-resolution obser- 
vations from pre-Swift GRB (e.g. Kulkarni et al. 1999; 
Barth et al. 2003) and the first year of Swift operation 
(e.g. Foley et al. 2005; Berger et al. 2005; Fynbo et al. 

2006) . Observers are now progressing toward the exami- 
nation of distributions of gas properties in the interstellar 
medium (ISM) of GRB host galaxies. 

Because of the apparent faintness of many GRB af- 
terglows and also instrument availability, the majority of 
GRB spectra are acquired with low-resolution spectrom- 
eters (FWHM> lA) and at moderate to poor SNR (< 
50 per resolution element). In these cases, abundance 
analysis must be performed on the equivalent width mea- 
surements of unresolved metal-line transitions. To date, 
the standard practice has been to perform a multi-ion 
single-component curve-of-growth (MISC-COG) analysis 
(e.g. Savaglio, Fall & Fiore 2003) or single-component pro- 
file fits to unresolved data (e.g. Savaglio & Fall 2004; Wat- 
son et al. 2005) which is fundamentally the same analysis. 
In this analysis, one derives an effective Doppler parame- 
ter bef f which principally characterizes the velocity width 
of strong, heavily saturated transitions. The column den- 
sity is then constrained with weaker, yet still potentially 
saturated transitions. In contrast, the H I column density 
can be determined from fits to the damping wings of the 
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Lya profile even with low-resolution data (e.g. Prochaska, 
Herbert-Fort & Wolfe 2005). 

The curve-of-growth technique dates far back (e.g. Unsold, 
Struve & Elvey 1930; Wilson 1939) and was primarily in- 
troduced to perform abundance analysis on lower spectral 
resolution data. The literature on the pitfalls of MISC- 
COG analysis is extensive (Nachman & Hobbs 1973; Crutcher 
1975; Spitzer & Jenkins 1975). These include systematic 
differences in the results for ions of differing ionization 
state (e.g. Mg I vs. Fe II; Spitzer & Jenkins 1975), 'hid- 
den' saturation due to significant variations in the Doppler 
parameter of overlapping components (Nachman & Hobbs 
1973), and the dangers of mixing refractory and non-refractory 
elements (Jenkins, Savage & Spitzer 1986). The Galac- 
tic ISM community is well aware of these issues; it has 
taken every effort to obtain high resolution spectroscopy 
and when limited to low resolution data, very cautiously 
interpret the observations. In short, few of the problems 
discussed in this paper are unique to GRB research, but 
we note that the discussion may translate to other abun- 
dance analyses with large, modern datasets of low resolu- 
tion spectroscopy (e.g. Savaglio et al. 2004; Turnshek et 
al. 2005). 

Jenkins (1986) examined the accuracy of a single-component 
COG analysis on multi-component line-profiles for a range 
of reasonable distributions of column density and Doppler 
parameter. He showed that the COG results typically give 
results accurate to within 20% (i.e. < 0.2 dex) of the cor- 
rect column density provided the N and b distributions of 
the individual components are well-behaved and so long 
as the peak optical depth tq < 5. Jenkins also empha- 
sized, however, that the results would likely significantly 
underestimate the column densities if the N oi b values 
were bi-modally distributed (see also Nachman & Hobbs 
1973). We will find that this is frequently the case for 
GRB line-profiles. 

In this paper, we will examine the main issues related 
to the ISM abundance analysis of GRB host galaxies as 
derived from low-resolution spectroscopy, i.e. equivalent 



1 



2 



On the Perils of Curve-of-Growth Analysis 



width analysis. We begin with an analysis of the GRB 051111 
sightline where Keck/HIRES observations permit one to 
directly test the MISC-COG technique (§ 2). The dis- 
couraging results are enlightening, namely the COG analy- 
sis systematically underestimates the column density mea- 
surements. We consider an additional example from the 
literature (GRB 020813) and raise similar concerns (§ 3). 
Finally, we offer a set of guidelines to consider when ana- 
lyzing low-resolution GRB spectra (§ 4). In a future paper 
(Prochaska et al. 2006), we will apply these guidelines 
to our datasets and previously published results with the 
intent of providing a uniform, database of measurements 
with realistic error estimates. 
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Fig. 1.— Fe II profiles from the ISM foreground to GRB 051111 
in the original data (LHS) and smoothed to lower resolution (RHS; 
FWHM f5! 5A) to mimic the line-profiles frequently observed. The 
velocity v = corresponds to z = 1.54948. The rest equivalent 
width values labeled on the Figure correspond to the full line-profile. 



2. GRB 051111: A TEST CASE 

To illustrate the key issues related to deriving abun- 
dances from equivalent width measurements, we will ex- 
amine a GRB sightline where the column densities are 
well-constrained, i.e., an afterglow with high quality, high 
resolution observations. The Keck observatory staff ac- 
quired HIRES (Vogt et al. 1994) observations of the bright 
afterglow of the z « 1.55 GRB 051111 (Hill et al. 2005) 



at high spectral resolution FWHM « 6 km s~^ and mod- 
erate signal-to-noise ratio (w 7 per 1.3km s~^ pixel) and 
released the data to the public^. The data reveal the de- 
tection of over 50 transitions (Prochaska et al. 2006) in- 
cluding a wealth of fine-structure levels (Prochaska, Chen 
& Bloom 2006). Our treatment begins with a single- 
component COG analysis of the Fe"*" ion and then pro- 
ceeds to a MISC-COG analysis of other ions in the ISM 
surrounding GRB 051111. 



Table 1 

EQUIVALENT WIDTH MEASUREMENTS FOR 
GRB 051111 
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> 13.64 


> 8.5 
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229.5 ± 4.1 
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213.2 ± 5.3 




> 4.7 



"J value for ions with excited states. 
''Rest equivalent width. 
'^Weighted-mean value. 

''Peak optical depth of the line profile assuming the column 
density from Column 6 and h = 7km s~^ . 



2.1. Single- component COG Analysis 

We will focus first on the Fe II transitions (Figure I) 
from the ground-state (J = 9/2) level and will adopt 
the Fe+ column density N(Fe+) = ioi5-27±o.oi ^^^^-2 
Prochaska, Chen & Bloom (2006). We have measured the 
total rest equivalent width of each Fe II transition and its 
associated statistical uncertainty (Table 1). At the reso- 
lution and signal-to-noise of these HIRES data, the rest 
equivalent width statistical uncertainty at 6OOOA in one 
pixel is a{W\) = 1.5mA. To better match the low res- 
olution observations, consider an analysis assuming total 
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uncertainty a(Wx) = 30mA. Figure 2 presents a single- 
component COG analysis of the Fe II transitions which is 
a two parameter fit to the data: the column density Nbest 
and the effective Doppler parameter b^f / • The upper panel 
shows the reduced equivalent widths W\ / A plotted against 
the /A product. The single-component model which min- 
imizes is overplotted on the observations. The model 
is driven to match the high W\ transitions because for 
fixed a{W\) these measurements have the highest relative 
SNR. This is a generic result for the analysis of unresolved 
line-profiles. The lower panel presents the Ax^ contours 
for the N, b^ff parameter space. The best-fit model (and 
its 2(7 error ellipse) underestimates the correct total Fe"*" 
column density by 0.6 dex, i.e. a factor of 4. 
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Fig. 2. — Single-component COG analysis of the Fe II transitions 
for GRB 051111 under the assumption that the equivalent width 
error cj{W\) = 30mA for each transition. The upper panel presents 
the Wx measurements and the model which minimizes x^- The 
lower panel shows Ax^ contours for the fee//, N parameter space. 
The best-fit column density (and 2cr error ellipse) severely underes- 



timate the true column density Ntru 



101 



Why does the single-component COG model systemat- 
ically underestimate the Fe^ column density? The princi- 
pal reason is that a single-component line-profile is only 
a good representation of the two weakest Fe II transi- 
tions (Figure 1). In order to match the observed equiv- 
alent width of the strong Fe II transitions which arc com- 
prised of many (> 10) 'clouds', a single-component COG 
analysis is driven to a large effective Doppler parameter, 
i.e. beff > 50kms^^. The only physical significance of 
this large 6e// value is that it describes the total veloc- 
ity width of the very strongest lines {Av ~ 200 km s~^). 



For this sightline, over 90% of the Fe+ ions are associ- 
ated with the component at u = Okms^^ and confined 
to Av w 20kms~^. Meanwhile, the gas responsible for 
> 80% of the equivalent width accounts for less than 10% 
of the total column density. It is this bi-modality in com- 
ponent column densities and the inclusion of highly satu- 
rated transitions which leads to such poor results for the 
COG analysis. 
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Fig. 3. — Offset between the best-fit column density Nf^^^t for a 
single-component COG analysis and the true column density Ntrue- 
In this analysis, we have adopted the Fe"*" optical depth profile from 
the unsaturated regions of the Fe II 2260, 2374, and 2586 transitions 
for GRB 051111. We then varied the normalization of the optical 
depth profile Tnrm, measured the W\ values for each of the transi- 
tions in Figure 1, and minimized for Afje^t and fee//- The solid 
curve shows the offsets assuming a(Wx) = 0.03mA for each transi- 
tion. The dotted curves show the results assuming cr{W\) is twice 
and four times smaller for the weakest transitions (Fe II 2249 and 
2260). Finally, the dashed lines show the results assuming cr(Wx) is 
2x and 4x larger for the weakest transitions. 



The best-fit Nbest value is sensitive to the error estimate 
a{W\) that we adopt, although only if we alter the rela- 
tive uncertainty. For example, if we set (t(W\) = 40mA for 
each transition then we will derive the same results as in 
Figure 2, albeit with a larger error ellipse. This is because, 
the relative of each line is unaltered. If we instead we 
increase a{W\) for only the strongest transitions, the best- 
fit model slowly approaches the correct value. Figure 3 
describes these trends. For this Figure, we have first de- 
rived the optical depth profile of the Fe"*" gas by combining 
the unsaturated observations of the Fe II 2260, 2374, and 
2586 profiles. We then varied the normalization of the op- 
tical depth profile (r„rm = 1 corresponds to the data as 
observed), measured the equivalent widths of the transi- 
tions, and minimized to derive Nbest and beff- The 
solid curve shows the results assuming a(Wx) = 30mA for 
each transition. At small Tnrm, the transitions are opti- 
cally thin and one recovers the correct result independent 
of the best-fit beff value. At Tnrm > 1, the Nbest value 
is systematically lower than Ntme- The dotted curves in 
the Figure show the results when we reduce the uncer- 
tainty in the Fe II 2249,2260 transitions by a factor of 2 
and 4 respectively. This leads to a significant improvement 
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for Tnrm < 5 where the hnes are at most modestly satu- 
rated. The dashed lines show the results where a{W\) is 
increased by 2 and 4 times for the weak transitions. The 
key conclusions to be drawn from Figures 2 and 3 is that 
the standard single-component COG analysis may lead to 
a severe underestimate of the gas column density. 

There are three signatures that the single-component 
COG analysis is providing an underestimate of the col- 
umn density. First, the effective Doppler parameter ex- 
ceeds 50km s~^ . Low-ion line-profiles in the Galactic ISM 
(e.g. Spitzer & Fitzpatrick 1995; Howk, Savage & Fabian 
1999), damped Lya systems (e.g. Prochaska & Wolfe 1996; 
Dessauges-Zavadsky et al. 2003), and GRB sightlines (e.g. 
Chen et al. 2005) all break into individual components 
with Doppler parameters h < 20kms^^. This is espe- 
cially true for the cloud which corresponds to the peak 
optical depth of the line. On its own, however, the deriva- 
tion of a large bejf value does not require that the mea- 
surement is erroneous Jenkins (1986). But, if one adopts 
beff 3> 20kms^^ then one will calculate column densi- 
ties for weaker transitions in the optically thin limit even 
if these transitions are actually saturated. A second sig- 
nature is that the reduced xt is large. This is likely to 
be a universal result for GRB observations, again because 
the single-component model is a poor representation of 
the data. The key concern here is that the error estimate 
that one derives from a A^^ analysis will generally by far 
too optimistic. Finally, the most important signature is 
when the reduced equivalent width of the weakest transi- 
tions are significantly under-predicted by the model. This 
is the most striking feature in Figure 2. We note that this 
effect is also evident, for example, in the Fe"*" measure- 
ments for GRB 000926, GRB 010222, and GRB 050505 as 
analyzed by Savaglio, Fall & Fiore (2003) and Berger et 
al. (2005). 
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Fig. 4. — Velocity profiles for a weak and strong transition arising 
in the ISM of the host galaxies of (a) GRB 050730, (b) GRB 050820, 
(c) GRB 051111, and (d) GRB 060418 (Prochaska et al. 2006). 
Note that the majority of column density is associated with a sin- 
gle, narrow component whereas the equivalent width of the strongest 
transitions is dominated by weaker 'clouds'. It is this bi-modal as- 
pect of the column density distribution which accentuates the fail- 
ings of the single component COG analysis. 



Before proceeding to the multi-ion, single-component 
COG analysis, we wish to comment on our results in light 
of the analysis by Jenkins (1986). As noted in the Intro- 
duction, Jenkins (1986) warned that the single-component 
COG analysis would significantly underestimate the col- 
umn density in cases where the optical depth distribution 
for the individual components is not smoothly distributed 
(i.e. bi-modal). Indeed, this is the case for the line- profiles 
for GRB 051111 (Figure 1). We emphasize that this is also 
the case for the ISM gas associated with GRB 050730 and 
GRB 050820, i.e. every case where moderate SNR, high- 
resolution echelle data exists (Chen et al. 2005; Prochaska 
et al. 2006). Again, the total optical depth of the line- 
profiles is dominated by a single, narrow component and 
the total equivalent width (for strong transitions) is dom- 
inated by low column density clouds. We expect that a 
bi-modal distribution of N values will be characteristic of 
many GRB sightlines'^: the high column density compo- 
nent may be related to the surrounding star-forming region 
and the lower column density gas may arise from the am- 
bient ISM or even halo of the galaxy. Figure 4 compares 
one weak and one strong transition from the ISM of the 
host galaxies of GRB 050730, GRB 050820, GRB 051111, 
and GRB 060418 (Prochaska et al. 2006). It is evident 
that the majority of column density is in a single narrow 
component while the equivalent width is dominated by 
multiple, weaker components. Of course, additional high- 
resolution observations of GRB afterglows will help reveal 
the typical N and b distributions of the line- profiles. 

2.2. Multi-ion Single- component COG Analysis 

Let us now consider a multi-ion single-component COG 
(MISC-COG) analysis for the GRB 051111 transitions. To 
date, this has been the standard approach for analyzing 
low-resolution GRB spectra. In this analysis, one min- 
imizes the global foi' ^ single b^ff value and unique 
Nbest values for each ion considered. Figure 5 shows the 
results for a number of ions observed toward GRB 051111. 
As before, we have assumed a-{W\) = 30mA for all of the 
transitions. Aside from the weakest Mg I and Fe II tran- 
sitions, we find the data is reasonably well represented by 
this model, primarily because each ion has an additional 
degree of freedom associated with it. In fact, one could eas- 
ily choose to ignore the 'discordant' Mg I and Fe II points 
on account of line-blending and/or an unidentified cosmic- 
ray (especially if the observations include only one or two 
examples). We emphasize, however, that these data points 
are severely discrepant because (i) the analysis includes 
at least one strong, heavily saturated transition (tq > 5) 
from each ion which has driven Ntest to a systematically 
low value; and (ii) the transitions are optically thin and 
therefore highlight that the Nbest value is too low. 

The lower panel of Figure 5 shows the offset between 
Nbest and Ntme for each of the ions. The effects for Si"*", 
Zn"*", and Mn"*" are expressed as upper limits (i.e. the off- 
set is at least as large as indicated) because the Ntme 
values derived from the HIRES observations are conserva- 
tive lower limits due to line saturation. It is evident that 
the abundances for all of the ions are systematically un- 

■^In passing, we also note that a bi-modal column density distribu- 
tion is also characteristic of many of the damped Ly« systems (e.g. 
Prochaska et al. 2001). 
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Fig. 5. — The multi-ion single-component curve-of-growth 
(MISC-COG) analysis for a series of ions observed in the ISM of 
GRB 051111. The upper panel shows the best-fit model and the 
observations. The lower panel shows the offset for the best-fit col- 
umn density and the actual value. The arrows indicate lower limits 
to the offsets because the Si II, Zn II, and Mn II transitions are all 
saturated at echelle resolution. 



N = 



Wx 
A 



/(A/A) 



(1) 



with A in Angstroms. To summarize, the analysis is driven 
by the strong Fe II transitions to a large b^f / value which 
leads to a systematic underestimate of the column density 
of all the ions. As worrisome, the underestimate varies 
from ion to ion and one cannot even calculate relative 
abundances to uncertainties of less than 0.3 dex. 

Not surprisingly, the results are improved if one elimi- 
nates the four strongest Fe II transitions from the analysis 
(Figure 6). We find, however, that we still derive a b^ff 
value that is too large and therefore significantly underesti- 
mate the column densities of several of the ions. If we only 
include those transitions which are unsaturated at echelle 



Fig. 6. — Same as Figure 5 except we have removed the strongest 
Fe II transitions from the analysis. The results are significantly 
improved due to the lower fee// value, yet the column densities of 
Mg I, Si II, Zn II, and Mn II are still significantly underestimated. 



derestimated with offsets of 0.15 dex to at least 0.7 dex. 
The case of Zn"*" is particularly problematic. The mea- 
surements consist of two Zn II transitions with /A values 
that differ by 0.3 dex. The observations are well fit by 
the MISC-COG analysis yet 7V(Zn+)^gj,j is at least 0.5 dex 3 Qj^g 020813- 
too low. The explanation is relatively simple: the effective 
Doppler is so large that these saturated profiles are treated 
as if they lie on the linear portion of the COG, i.e.. 



resolution^ (i.e. consider only those lines which are truly on 
the linear portion of the COG), one finally achieves reason- 
able results. With low-resolution spectroscopy, however, 
one does not have this luxury. We are particularly wor- 
ried that low-resolution spectroscopy is generally limited 
to the analysis of transitions with \og{Wx/ X) > —4. 



A HIGH SNR, LOW-RESOLUTION EXAMPLE 
FROM THE LITERATURE 



Let us now consider a low-resolution spectroscopic ex- 
ample from the literature, GRB 020813. Barth et al. 
(2003) obtained a high signal-to-noise, low-resolution spec- 
trum of this afterglow with the LRIS spectrometer (Oke 
et al. 1995) in spectropolarimetery mode. Barth et al. 
(2003) reported equivalent width measurements for many 
transitions associated with the ISM surrounding GRB 020813. 
We present in Figure 7 a single-component COG analysis 
of the Fe II transitions detected in their data. Similar 
to GRB 051111 (Figure 2), the strong transitions drive 
the best-fit model to a large beff value. In contrast to 
GRB 051111, however, the weak Fell 2249 and Fell 2260 
transitions are well- modeled by this fit. The error el- 
lipse for this model indicates a very precise measurement, 



A^(Fe+) = 10 



115.52±0.03 



cm , in good agreement with the 



line-profile analysis of Savaglio & Fall (2004). The corre- 
spondence between the two analyses does not lend addi- 
tional confidence; for unresolved features, a single-component 

■^This does not include Si+ or Zn+ whose lines have tq > 5. 
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Fig. 7. — Single-component COG analysis of the Fe II transitions 
for GRB 020813. The upper panel presents the W\ measurements 
and the model which minimizes \^ . The lower panel shows Ax^ 
contours for the be//i ^ parameter space (note the tiny dynamic 
range). 
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Fig. 8.— Model profiles for the Fe II transitions of GRB 020813. 
The LHS shows the single-component profile corresponding to the 
best-fit model for Figure 8. Note that the two Fe II transitions are 
optically thin even though they exhibit relatively large Wx values. 
The RHS shows a toy model of the line-profile consisting of two 
strong, narrow components and a series of weak features. It presents 
a good match to the observed equivalent widths yet has 100% higher 
column density than the single-component model. 



line-profile analysis is identical to the COG method. 

Because of the good visual agreement of the model (for- 
mally the value is very poor) and the very small error 
ellipse, one is led to believe that the column density has 
been precisely determined. This is an incorrect conclusion. 
First, consider the best-fit value Attest = 10^^'^^ cm^^. 
Adopting 6e// = 65 km s~^, the line-profile of the Fe II 2249 
and 2260 transitions are optically thin, i.e. the transitions 
are predicted to lie on the linear portion of the COG (Fig- 
ure 7). In this case, one can calculate A^(Fe"'") from Equa- 
tion 1 and find values of lO^'^ "'^ cm~^ and 10^^-^'' cm^^ for 
Fe II 2249 and 2260 respectively. While this result appears 
to be self-consistent, consider the implication: the rest 
equivalent widths of these two transitions are very large, 
in fact 2.5 X larger than the partially saturated Fe II tran- 
sitions for GRB 051111 (Figure 1). To account for these 
large equivalent widths and the best-fit column density, 
one must hypothesize a line-profile consisting of a series of 
non-overlapping components that have large column den- 
sity, but not too large [N ~ 10"'^^ cm~^). While this cannot 
be ruled out by the observations, we consider this scenario 
to be highly unlikely in light of the observed line-profiles 
of GRB sightlines at high-resolution (e.g. Figure 1). 

To emphasize this point, consider Figure 8 where we 
compare the single-component scenario with an alternate 
model. The profiles on the left-hand side show several 
of the observed Fe II transitions predicted by the best-fit 



model. The right-hand side shows a toy model consist- 
ing of two strong, narrow components (each with logiV = 
15.4 and 6 = 5 km s^^) and a series of weak components 
included to reproduce the large equivalent width of the 
stronger Fe II transitions. This toy model provides a rea- 
sonable match to the observed equivalent widths (in fact, 
a better fit to the strong transitions) and implies a total 
Fe"*" column density that is two times larger (0.3 dex or 
lOtr) than the single-component COG result. 

We do admit, however, that it may be difficult to con- 
struct an optical depth profile for GRB 020813 where the 
true Fe+ column density is 10 x larger than the single- 
component COG value. And, we note that the low SNR, 
high-resolution VLT/UVES observations of GRB 020813 
(Fiore et al. 2005) do suggest that the optical depth pro- 
file is not dominated by only a single strong component. 
Nevertheless, we consider it impossible that the best-fit 
COG value is not underestimating the N(Fe'^) value. Fur- 
thermore, we are certain that the statistical error estimate 
grossly underestimates the systematic uncertainty associ- 
ated with one's assumptions about the optical depth pro- 
file. 

Continuing with the treatment of GRB 020813, Figure 9 
presents the MISC-COG analysis for a series of the ions 
detected in the Barth et al. (2003) spectrum. Again, 
the model is visually a good fit to the observations and 
we derive column densities that are in good agreement 
with Savaglio & Fall (2004). The issues that exist for Fe+ 
are accentuated for the Si"*" and Zn+ results. Here, the 
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Fig. 9.— MISC-COG analysis of the ions in the ISM of 
GRB 020813. Because of the large be// values demanded by the 
heavily saturated Fe 11 transitions, the other transitions are treated 
as if they lie on the linear portion of the curve-of-growth. This leads 
to unrealistically low central values for the best-fit column densities 
and overly optimistic error estimates. 



rest equivalent widths are Wx = 0.7A and 0.6A for the 
Si II 1808 and Zn II 2026 profiles, i.e. very large. Yet, 
again, the best-fit column densities correspond to the op- 
tically thin hmit because of the large 6e// value. It is 
nearly impossible to construct optically thin Si II 1808 
and Zn II 2026 line-profiles if one is restricted to physi- 
cally realistic individual clouds with 6 < 10 km s~^. Fur- 
thermore, there is no observational constraint which sets 
an upper limit to the Si'^ and Zn^ column densities. 

These are generic results for any MISC-COG analysis 
that gives a very large 6e// value: the technique models 
transitions with very large equivalent width measurements 
as being optically thin and adopts the lowest column den- 
sity conceivable. The bottom line is that one is very likely 
to underestimate the column densities of all of these ions. 
And, the error estimates should reflect a one-sided distri- 
bution allowing only significantly larger values, i.e. lower 
limits. 

4. CONSERVATIVE SUGGESTIONS FOR EQUIVALENT 
WIDTH ANALYSIS 

Given the cautionary (pessimistic) tone of the previous 
sections, one may draw the conclusion that no progress 
on chemical abundances can be made with low-resolution 
GRB spectroscopy. To push forward in a productive yet 
cautious manner, we offer the following suggestions. 

1. Assume that every line-profile is saturated and ex- 
amine all possible constraints to demonstrate otherwise. 
Absent additional constraints, the equivalent width mea- 
surements establish only lower limits to the column densi- 
ties. 

2. Consider independent analyses of the heavily satu- 
rated and the weakest lines. The former may give esti- 
mates for the abundances of the ambient ISM along the 
sightline while the former could be more relevant to the 
gas near the star-forming region. Unfortunately, one can- 
not separate the H I gas along the sightline into these 



two phases so the goals of such analysis would be to (i) 
prevent the strong strong lines from biasing the total col- 
umn density to low values and (ii) attempt to measure 
relative abundances in the two 'phases'. We reemphasize 
that by considering one COG analysis for all the lines, the 
strong lines will drive the COG analysis to large 5e// val- 
ues which forces the column densities of weaker (yet likely 
saturated) transitions to their minimum value. The end 
result is a systematic underestimate of all column densi- 
ties. A corollary to this suggestion is that the results for 
a single-component COG analysis with 6e// > 20 km s~^ 
are highly suspect. 



Table 2 

WEAK LINES FOR UPPER LIMITS 



Transition 


A 


log/ 


TO = 2 

( cm ^) 


wl 


[M/H]': 




(A) 




(mA) 




SII 1250 


1250.5840 


-2.2634 


15.29 


82 


-1.91 


OI 1355 


1355.5977 


-5.9066 


18.90 


89 


0.16 


BII 1362 


1362.4611 


-0.0057 


13.00 


90 


0.21 


Nill 1467b 


1467.7560 


-2.0044 


14.96 


97 


-1.29 


PII 1532 


1532.5330 


-2.1186 


15.06 


101 


-0.47 


Gell 1602 


1602.4863 


-0.8428 


13.76 


106 


0.13 


Fell 1611 


1611.2004 


-2.8665 


15.79 


106 


-1.71 


Nill 1703 


1703.4050 


-2.2218 


15.12 


112 


-1.13 


Mgl 1747 


1747.7937 


-2.0419 


14.93 


115 


-2.65 


Fell 1901 


1901.7729 


-3.9961 


16.84 


125 


-0.66 


CrII 2066 


2066.1609 


-1.2882 


14.10 


136 


-1.57 


Fell 2249 


2249.8767 


-2.7397 


15.51 


148 


-1.99 


Sill 2335 


2335.1230 


-5.3716 


18.13 


154 


0.57 


Till 3230 


3230.1311 


-1.1630 


13.78 


213 


-1.16 



"Column density corresponding to a line-profile with peak optical depth 

TO = 2 and Doppler parameter b = 10 km s~^. 
''Rest equivalent width for the column density in column 4. 
■^Metallicity assuming the column density in column 4 and logAfjji = 22. 

3. High precision measurements, absolute or relative, 
will be very difficult to achieve. Meaningful constraints 
that require better than 0.5 dex precision (e.g. nucleosyn- 
thetic patterns) may not be attainable. One may, however, 
be able to distinguish between very high and low dust-to- 
gas ratios. 

4. Aim for the highest spectral resolution at the cost of 
high signal-to-noise. The only convincing argument that 
a line is unsaturated is to reconstruct the optical depth 
profile. 

5. Place upper limits on the metallicity using very weak 
transitions and/or trace elements. Table 2 provides a list 
of weak transitions for abundant elements (e.g. Fe, Si, O) 
and stronger transitions of trace elements (e.g. Ge, B). 
The Table also lists the column density corresponding to 
a line-profile which is modestly saturated (tq — 2) for a 
Doppler parameter 6 = 8 km s^^. If one can establish an 
upper limit (ideally 3cr including continuum uncertainty) 
to the equivalent width which is below the corresponding 
W\ value (Column 5) then one places an upper limit to 
the metallicity listed in Column 6, assuming solar relative 
abundances and ignoring differential depletion. In general, 
this will require observations with rest equivalent width 
error a{W\) < 30mA. 
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5. CONCLUDING REMARKS 

Before concluding, we wish to briefly comment on the 
abundance measurements for the ISM of GRB presented 
in the Uterature. Table 3 summarizes the results to date. 
The last two columns give the minimum and maximum 
metallieity values allowed by the data as reported in the 
literature^. One notes that solar (and even super-solar) 
abundances are allowed in nearly every case. It is evident 
from Table 3 that the principal conclusion of Savaglio, Fall 
& Fiore (2003) - GRB sightlines have extreme metal col- 
umn densities may be a remarkable understatement. We 
also caution that the common conclusion that GRB occur 
in a metal-poor ISM ([M/H] < — 1) may need serious re- 
vision. 

It is unfortunate that there are few transitions for mild 
or non-refractory elements that are weak enough to place 
meaningful upper limits to the abundances; the Zn and Si 
abundances, for example, are compromised by their rel- 
atively strong transitions. We expect, however, that an 
analysis of the lines listed in Table 2 may improve the 
constraints (Prochaska et al. 2006) and provide more re- 
alistic metallieity measurements for the gas within GRB 
host galaxies. 
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Table 3 

SUMMARY TABLE OF ABUNDANCE MEASUREMENTS IN LOW-RESOLUTION GRB SPECTRA 



GRB 




z 


log AThi 


[M/H]», 


Rcf 


Instr 


Line" , 


W^a(A) 


[M/H]^„ 




990123 


1, 


.600 


22 


~'J-'-0.05 


1,2 


Keck/LRIS 


Zn II 2026 


0.485 ± 0.038 


-0.95 




000926 


2, 


.038 


21.3 ±0.25 




2,3,4 


Keck/ESI 


Zn II 2026 


0.90 ±0.1 


-0.25 




010222 


1, 


.477 


22 




5,2 


Keck/ESI 


Zn II 2026 


0.91 ±0.1 


-1.0 




011211 


2 


.142 


20.4 ±0.2 


-0.91°:! 


6,2 


VLT/FORS2 


Zn II 2026 


0.91 ±0.1 


-1.0 




020813 


1, 


.255 


22 


1 1+0.06 

-0.06 


7,2 


Kcck/LRIS 


Zn II 2026 


0.61 ±0.1 


-1.2 




030323 


3, 


.372 


21.9 ±0.07 


1 O+0.20 
-'-••^-0.20 


8 


VLT/FORS2 


Zn II 2026 


0.82 ± 0.09 


-0.9 


±0.5 


050401 


2, 


.899 


22.5 ± 0.3 


-o.8±r4 


9 


VLT/FORS2 


Zn II 2062 




-1.0 




050505 


4, 


.275 


22.05 ± 0.1 


-1.2 


10 


Kcck/LRIS 


Fe II 1611 


0.54 ± 0.05 


-1.0 




050904 


6, 


.296 


21.3 ±0.2 




11 


Subaru/FOCAS 


S II 1253 


0.5 ±0.15 


-1.4 




060206 


4, 


.048 


20.85 ±0.1 


-0.84+°i 


12 


WHT/ISIS 


S II 1250 




-1.0 


-0.5 



"Metallicity reported in the literature. For those sightlines without ATjji measurements, we adopt a reasonably conservative value of 10^^ cm ^. 
Note that the errors reported reflect only the uncertainy in the reported metal abundance to highUght the overly optimistic results. 

"Most constraining transition for a mctallicity measurement, generally derived from a weak transition. Note that we have not corrected for Mg I 
blending with the Zn II 2026 transition. This could reduce the mctallicity limit by as much as 0.2 dex. 
'^Minimum metallicity derived from assuming the transition in Column 8 is optically thin. 

"^Maximum metallicity based on the non-detection of a transition or the accurate measurement of an unsaturated transition. 
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